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ABSTRACT 

We examine the low angular momentum flow model for Sgr A* using two-dimensional 
hydrodynamical calculations based on the parameters of the specific angular momen- 
tum and total energy estimated in the recent analysis of stellar wind of nearby stars 
around Sgr A*. The accretion flow with the plausible parameters is non-stationary 
and an irregularly oscillating shock is formed in the inner region of a few tens to a 
hundred and sixty Schwarzschild radii. Due to the oscillating shock, the luminosity 
and the mass-outflow rate are modulated by several per cent to a factor of 5 and a 
factor of 2-7, respectively, on time-scales of an hour to ten days. The flows are highly 
advected and the radiative efficiency of the accreting matter into radiation is very 
low, ~ 10 _5 -10 -3 , and the input accretion rate Mi npu t = 4.0 x 10~ 6 M Q yr -1 results 
in the observed luminosities ~ 10 36 erg s _1 of Sgr A* if a two-temperature model 
and the synchrotron emission are taken into account. The mass-outflow rate of the 
gas originating in the post-shock region increases with the increasing input specific 
angular momentum and ranges from a few to 99 per cent of the input accreting mat- 
ter, depending on the input angular momentum. The oscillating shock is necessarily 
triggered if the specific angular momentum and the specific energy belong to or are 
located just nearby in the range of parameters responsible for a stationary shock in 
rotating inviscid and adiabatic accretion flow. The time variability may be relevant to 
the flare activity of Sgr A*. 

Key words: accretion, accretion discs - black hole physics - hydrodynamics - radi- 
ation mechanism: thermal - shock waves - Galaxy: centre. 



1 INTRODUCTION 

Disc accretion is an essential process for such phenomena 
as energetic X-ray sources, acti ve galactic nuclei and p roto- 
stars. Since the early work s by iPringle fc Reesl (|l972l ^ and 
IShakura fe Sunvaevl l | 1973ft . a great number of papers have 
been devoted to studies of the disc accretion onto gravitating 
objects. When the accretion rate is not too high, the accre- 
tion disc luminosity is directly in proportion to the accretion 
rate and can be successfully described by the standard thin 
disc model (S-S model). Sgr A* in our Galactic centre has 
been extensively studied in the category of accretion pro- 
cesses because it is a supermassive black hole candidate in 
our Galaxy and has unique observational features incom- 
patible with the S-S model. One of the features of Sgr A* 
is that the observed luminosity is five orders of magnitude 
lower than that predicted by the S-S model. Moreover, the 
spectrum of Sgr A* differs from the multi-temperature black 
body spectrum obtained from the S-S model. 
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Since the observational features of Sgr A* can not be ex- 
plained by the S-S model, two types of theoretical models, 
such as the spherical Bo ndi accretion model without any 
net angular momentum (|Bondil [l952) and the advection- 
dominated accretion flow model (ADA F) with high angular 
momentum l|Naravan fc Yi|[l994l . 1 19951 ) , have been proposed 
(Narayan & McClintock 2008, Yuan 2011 for review). Both 
the Bondi model and the ADAF model result in highly ad- 
vected flows and the radiative efficiency is so low as to be 
compatible with the observations. The spherically symmet- 
ric Bondi model gave the first detail ed model of Sgr A* 
|Melialll992l : lMelia. Liu fc Cokerll200ll ). This model applies 
when the circularization radius is smaller than the inner- 
most stable circular orbit. The Bondi model is frequently 
used to estimate the expected accretion rate near the Bondi 
radius far from the gravitating object. However, the simple 
Bondi model can not be extended close to the black hole to 
satisfy the details of the activity of Sgr A*. In contrast with 
the Bondi model, the ADAF model applies when the angular 
momentum at the outer boundary is roughly Keplerian, and 
a large fraction of the viscously generated heat is advected 
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with the accreting gas and only a small fraction of the energy 
is radiated. Unlike the cold standard disc, this model results 
in a high temperature, rapid accretion, sub-Keplerian disc 
and inefficient radiation in the inner region of the disc. The 
ADAF models were sho wn to be generally successful and 
more advanced models l|Yuan. Quataert fc Naravanl 120031 . 
12004 ) . taking into account the parametric description of the 
outflow and jet, explain well the observations. The impor- 
tant key to these models for Sgr A* is the amount of angu- 
lar momentum in the accretion flow. However, at present, 
we have no clear evidence for the angular momentum from 
observations. 

The low angular momentum model belongs to an inter- 
mediate case between the Bon di model and the ADAF model 
and was applied to Sgr A* (IMoscibrodzka. Das fc Czernvl 
120061 ; ICzernv fc Moscibrodzkal I2008T ). Assuming that the 
Wolf-Rayet star IRS 13 E3 is the dominant source of the 
matter accreting onto Sgr A* and assuming the wind tem- 
perature T w i n( j = 1.0 or 0.5 keV, they estimated the net 
angular momentum A of 1.68-2.16 and the Bernoulli con- 
stant s of 3.97 x 10" 6 -1.98 x 10" 6 , where the mass M, the 
speed of light c and the Schwarzschild radius R g = 2GM/c 2 
are used as the units of mass, velocity and distance. With 
these flow parameters for A and e, they showed analytically 
that there is no continuous flow solution which attains to 
the event horizon, and the resulting flow would be non- 
stationary, but that, for the case of the angular momentum 
A= 1.55 lower than the best estimates for Sgr A*, there ex- 
ists a standard stationary shock solution. Subsequently, they 
speculated that, for the non-stationary flow, the inflowing 
matter accumulates and forms a ring or torus and that the 
ring in the inner region is unstable and may explain the vari- 
ability of Sgr A*. Motivated by their suggestion and results, 
we examine the low angular momentum flow model for Sgr 
A* using 2D hydrodynamical calculations and discuss the 
results on the activity of Sgr A*. 



2 MODELLING A LOW ANGULAR 
MOMENTUM FLOW 

We consider a supermassive black hole with mass M = 
4 x 10 6 M Q for Sgr A*. For the accretion flow around Sgr 
A*, we use the flow parameters A and e estimated by 
IMoscibrodzka. Das fc Czernvl l|2006h and the mass accretion 
rate M = 4.0 x 10~ 6 M Q yr _1 in all models. ID flow solution 
of a thin, rotating, inviscid and adiabatic accretion flow near 
the equatorial plane generally depends on the flow model 
especially on its thickness, for example, the vertical equilib- 
rium model or the constant height model. Assuming the ver- 
tical equilibrium model here, we solve the Bernoulli equation 
with the constant angular momentum A, find the outer and 
inner sonic points and get the Mach number versus radius 
relation, the sound speed v s , the thickness h of the accretion 
flow, the radial velocity v an d the Temperature T at a given 
radius r l|Chakrabartilll989l ) . It should be noticed that the 
configuration of the Mach nu mber versus radius does not co- 
incide completely with that in IMoscibrodzka. Das fc Czernvl 
(2006) because of the different flow model used. The Mach 
number versus radius relations for A = 1.68 and 2.16 are 
similar to theirs, that is, the flow which passes through the 
outer sonic point never continue to the event horizon. The 



flow for A = 1.55 also shows a flow topology similar to that 
for A = 1.68, unlike the case of lMoscibrodzka. Das fc Czernvl 
(2006). The first purpose in this study is to examine how the 
flow with these specific angular momenta actually behaves. 
However, in addition to the above parameters, we examined 
another case of A = 1.35, for which the standing shock is 
predicted to exist in the inner region. 

Figs 1 and 2 show the Mach number versus radius in 
ID adiabatic flow obtained for the cases of A = 1.68 and e 
= 3.97 x 10~ 6 (model A) and A = 1.35 and e = 1.98 x 10~ 6 
(model D), respectively. In Fig. 1, the particle which passes 
through the outer sonic point A falls down supersonically 
inward but never attains the event horizon since it makes a 
closed loop of the Mach number curve. On the other hand, 
in Fig. 2, the particle falls supersonically along the solid line 
AB, jumps at the shock position of R s ~ 20-R g from the 
point B to C in a subsonic state and tends supersonically 
toward the event horizon along the dashed line CD, where 
the Mach number curve (dashed line) which passes through 
the inner sonic point intersects at the point C with the other 
Mach number curve AC (dotted line) which passes through 
the outer sonic point and satisfies the Rankine-Hugoniot re- 
lation between the solid and the dashed lines. In Fig. 2, there 
exist two shock positions (R s ~ 17 and 19). It is known that 
usually the inner sh ock of the two sh ocks is unstable but the 
outer one is stable (|Nakavamall 19941 ). 

Thus far, the model parameters are listed in Table 1, 
where 7 is the adiabatic index. First we examine models 
A-D using 2D radiation hydrodynamical calculations which 
include only the free-free transitions in the cooling and heat- 
ing rate. It is well known that the existence of the sub- 
millimetre bump in Sgr A* spectra is certainly caused by the 
synchrotron emission. Secondly, following the results of mod- 
els A and D, we examine another models of E and F which 
take account of the synchrotron cooling in the energy equa- 
tion under an assumption that the magnetic energy density 
is 10 per cent of the thermal energy density. In models E 
and F, we adopt a two-temperature model where the ratio 
Te/Ti of the electron temperature T e to the ion temperature 
Ti is initially assumed to be a constant 1) and neglect 
the radiation transport because the flow is sufficiently op- 
tically thin and the flux-limitted diffusion approximation of 
radiation is less accurate for so optically thin medium. 



3 MODEL EQUATIONS 

The set of relevant equations consists of six partial differ- 
ential equations for density, momentum, and thermal and 
radiation energy. These equations include the heating and 
cooling of gas and radiation transport. The radiation trans- 
port is treated in the grey, flux-limit ed diffusion approxima- 
tion dLevermore fc Pomra ning 198l]). We use spherical polar 
coordinates (r,(,ip), where r is the radial distance, £ is the 
polar angle measured from the equatorial plane and y> is the 
azimuthal angle. The flow is assumed to be symmetrical with 
respect to the Z-axis {d/d<p = 0) and the equatorial plane. 
In this coordinate system, we have the ba sic equations in 
the following conservative form (lKlevlll989l ): 



dp 
dt 



+ div(pv) — 0, 



(1) 
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Table 1. Model parameters for the accreting matter onto Sgr A*, where in models E and F a two-temperature model is adopted and 8 
is the initial ratio T e /T; of electron temperature T e to ion temperature T;. 



model 


A 


e 




7 


M (M Q yr" 1 ) 


5 


A 


1.68 


3.97 xlCT 


ii 


1.6 


4.o xicr 6 




B 


2.16 


1.98 XHT 


-6 


1.6 


4.0 xlO -6 




C 


1.55 


1.98 xlO" 


-6 


1.6 


4.0 xl0~ 6 




D 


1.35 


1.98 Xl(T 


-6 


1.6 


4.0 xl0~ 6 




E 


1.68 


3.97 xlCT 


-6 


1.6 


4.0 xl0~ 6 


3/7 


F 


1.35 


1.98 xlCT 


-6 


1.6 


4.0 xl0~ 6 


1/9 



d(pv) 
dt 



+ div(pvv) = p 



+ 



GM 



(r-P g 



dp 
dr 



d(prw) 2 9p 

— — h dw(prwv) = -pu^tanC - ^- + fc, 



+ /r,(2) 



(3) 



d(pr cos(v v ) 

di 



+ div(prcos£v ip v) = 0, 



(4) 
(5) 



+ div(pev) = — p divv — A, 
and 

^ + divFo + div(vE + v • P ) = A - p ( K + a K ■ F , (6) 

at c 

where p is the density, v = (v, w,Vcp) are the three velocity 
components, G is the gravitational constant, p is the gas 
pressure, e is the specific internal energy of the gas, Po is 
the radiation energy density per unit volume, and Po is the 
radiati on stress tensor. We adopt the pseudo-Newtonian po- 
tential (|Paczvriskv fc Wiitalll980t ) in equation (2). The force 
density / R = ( f r , fc ) exerted by the radiation field is given 
by 



In = P -Po, 



(7) 



where k and a denote the absorption and scattering coeffi- 
cients and Fo is the radiative flux in the comoving frame. 
The quantity A describes the cooling and heating of the gas, 



A = pcn(S* - E 



(8) 



where S* is the source function. For this source function, 
we assume local thermal equilibrium S* = aT 4 , where T is 
the gas temperature and a is the radiation constant. For the 
equation of state, the gas pressure is given by the ideal gas 
law, p — RgpT ' I p, where p is the mean molecular weight and 
Rg is the gas constant. The temperature T is proportional to 
the specific internal energy, e, and satisfies the relation p — 
(7 — l)pe = RGpT/p. To close the system of equations, we 
use the flux-limited diffusion approximation for the radiative 
flux: 



Fo = - 



and 



A c 



p(K- 



-grad Eo, 



Pq = Eo ■ TEdd, 



(9) 



(10) 



where Ao and TEdd are the flux-limiter and the Eddington 
Tensor, respec tively , for w hich we use the approximate for- 
mula given in iKlevI (|l989h . The formula fulfills the correct 
limiting conditions Ao — > 1/3 in the optically thick diffusion 



limit, and | Fo \—¥ cEo as Ao —5- for the optically thin 
streaming limit. 

In the two-temperature model with Eo=0 for models E 
and F, the energy equation (5) is replaced by the following 
energy equations for ion and electron. 



+ div(peiv) = -pi divv - q ie , 

and 
dpe c 



dt 



div(pe c v) = —p c divv - pc^S, + qi c — q syn , 



(11) 



(12) 



where e; and e e are the specific internal energy of the ion 
and electron, pi and p e are the ion and electron gas pressure, 
qi c is the energy transfer rate from the ion to electron by 
Coulomb collisions, fjr Bvn is the cooling rate by synchrotron 
radiation (Na ravan fc Yilll995T ). and p = p\ + p c . 



4 NUMERICAL METHODS 

The set of partial differential equations (l)-(6) is numeri- 
cally solved by a finite-difference method under adequate ini- 
tial and boundary conditions. The numerical schemes used 
are ba sicall y the same as those previou s ly des cribed in IKlevI 
l|l989h and fokuda. Fuiita fc Sakashital (|l997l 1. These meth- 
ods are based on an explicit-implicit finite difference scheme. 
The inner boundary radius R[ n is taken to be 27? g . In our 
models, the outer shock will be formed in the region of 
r ^ 200 R s if it exists. Therefore, the outer boundary ra- 
dius Pout is set to be 200P g so that the predicted shock 
position never exceed the outer boundary. The computa- 
tional domain is divided into N r x grid cells, where A^ r 
grid points (=100) in the radial direction are spaced loga- 
rithmically as Ar/r — 0.054 and A^ grid points (=100) in 
the angular direction are equally spaced. 



4.1 Initial flow variables 

The initial flow variables in the 2D time-dependent calcula- 
tion for models A-D are given from the ID adiabatic flow 
mentioned in Section 2. To obtain the radiation energy den- 
sity Po, we specify the ratio /3 of the gas pressure to the total 
pressure and the Eddington factor f^dd- Then, the gas tem- 
perature T and the radiation energy density Po are given 
by 



T = /3 Us 2 /(Pg7), 
and 

Eo = (1-/3KV(/ E dd7), 



(13) 



(14) 
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Figure 1. The topology of Mach 
number versus radius of the 
transonic flow for the parame- 
ters of model A, where the solid 
and the dashed lines intersect 
at the outermost and innermost 
sonic points indicated by crosses, 
respectively. The particle which 
passes through the outer sonic 
point A falls down supersonically 
but never attains the event hori- 
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Figure 2. Same as Fig.l but 
for model D. The particle which 
passes through the outer sonic 
point A falls supersonically 
along the line AB, jumps at 
the shock position R B ~ 20R g 
from the point B to the point 
C in a subsonic state and again 
tends towards the event horizon 
supersonically. 



Table 2 shows the flow variables of the radial velocity « out , 
the Mach number M ac h, the density p out , the gas tempera- 
ture Tout, the radiation energy density (-Eo)outj the relative 
thickness h/r of the flow at the outer boundary flout, the ra- 
dius -Raonic of the outer sonic point and the gas temperature 
Tb at the Bondi-Hoyle-Lyttleton radius. 

On the other hand, as for the initial flow of models 
E and F, we use the final numerical results of models A 
and D, respectively and assume initially that T e /T[ = S and 
T — (T c +Ti)/2 to maintain the pressure balance throughout 
the flow, where T is the single temperature obtained from 
models A and D. 



4.2 Initial and boundary conditions 

At the outer boundary of the accretion flow, a continuous 
inflow with the flow variables in Table 2 is always given. The 
initial atmosphere above the accreting matter in the com- 
putational domain is given as a radially hydrostatic equilib- 
rium state with zero azimuthal velocities everywhere. Physi- 
cal variables at the inner boundary, except for the velocities, 
are given by extrapolation of the variables near the bound- 
ary. However, we impose limit conditions that the radial 
velocities are given by a free-fall velocity and the angular 
velocities are zero. On the rotational axis and the equatorial 
plane, the meridional tangential velocity w is zero and all 
scalar variables must be symmetric relative to these axes. 
As for the outer boundary region above the outer accreting 
matter, we use free-floating conditions and allow for the out- 
flow of matter, whereas here any inflow is prohibited. With 
these initial and boundary conditions, we perform time in- 
tegration of equations (l)-(6) until a quasi-steady solution 
is obtained. 



5 NUMERICAL RESULTS 

To get a steady state of flow, we need evolutionary times 
longer than the following dynamical time tdyn- 



4 x 10 4 / 



A 



V 200i? g 



(15) 



where Q is the angular velocity. 

This time is a measure of the needed computational 
time. Actually, because unsteady shock phenomena oc- 
curred, we needed more times, of 5-7 x 10 4 Rg/c, to examine 
the time variabilities, where R g /c ~ 40 s. The total lumi- 
nosity L is also a good measure to check whether a steady 
state of the accretion flow is attained. The luminosity L is 
given by 

fTT/2 

L = 4xRo^ I F r {R out ,C,)cos(d(. (16) 



Concerning /?, we specify its value so that the gas pres- 
sure is dominant and the temperature Tb at the Bondi- 
Hoyle-Lyttleton radius -Rbhl (~ 7 x 10 4 i? g ) is comparable 
to the assumed wind temperature T w , where Rbhl is given 
by the wind vel ocity 1000 km s -1 and the wind tem perature 
T w = 1.0 keV (|Moscibrodzka. Das fc Czenry||2006ft . Conse- 
quently, /3 = 0.9 and /Edd = 1 are adopted in all models. 



5.1 Time variations of the luminosity and the 
mass— flow rate 

First, we show the results of model A as a typical case. After 
a time of ~ 10 3 7? g /c, the total luminosity L through the 
outer boundary nearly attains a steady value of 5.4 x 10 38 
erg s _1 but is modulated only by a few to ten per cent of the 
total luminosity. Following the time evolution of the flow, a 
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Table 2. Flow variables at outer boundary R ou t, radius -R aon ic of the outer sonic point, gas temperature T B at the outer sonic point and 
gas temperature Tb at the Bondi-Hoyle-Lyttleton radius -Rbhl, where po = 10~ 8 g cm -3 . The finally obtained flow variables of models 
A and D are used as the initial flow variables of models E and F, respectively, but with ion temperature (T;) ou t = 2T ou t/(l + S) and 
electron tempcratute (T c ) out = 2<5T out /(l + 5). 



Model 


Vout/c 




M ach 


Pout 




Tout 


(^o)out 




(h/r) out 


-^sonic 


T s 












(g cm~ 


3 ) 


(K) 


(POC 2 ) 






(Kg) 


(K) 


(K) 


A 


4.75x10" 


-2 


1.91 


5.86x10 


-19 


2.45 xlO 9 


8.86x10- 


15 


0.45 


8.359x10 s 


9.5X10 7 


9.6x10 s 


B 


4.88xirr 


-2 


2.02 


5.90x10 


-19 


2.30 xlO 9 


8.35x10- 


15 


0.43 


1.677X10 4 


4.7X10 7 


l.OxlO 7 


C 


4.95x10" 


-2 


2.06 


5.83x10 


-19 


2.28 x10 s 


8.20x10- 


15 


0.43 


1.681X10 4 


4.7X10 7 


l.OxlO 7 


D 


4.97x10- 


-2 


2.07 


5.81x10 


-19 


2.28 xlO 9 


8.16x10- 


15 


0.45 


1.682X10 4 


4.7X10 7 


l.OxlO 7 



shock wave is formed in the inner region of r < 30i? g and 
it moves gradually outward. After a time of ~ 10 4 Rg/c, the 
shock wave begins to oscillate irregularly between the radii 
of 80-160 R g with time-scales of hours to days. The mass- 
outflow rate M ou t at the outer boundary and the mass- inflow 
rate Mi n at the inner edge are modulated intermittently by a 
factor around 3. M ou t plus M n is nearly equal to the input 
accretion rate Mi npu t but it is remarkable that the mass- 
outflow rate is twice as large as the mass-inflow rate at the 
inner edge, that is, a large amount of the input accreting 
matter is ejected as a wind. 

Fig. 3 shows the time variations of the total luminos- 
ity L (erg s" 1 ), the mass-outflow rate M out (g s _1 ), the 
mass-inflow rate M n and the shock position on the equa- 
torial plane for model A, where the input accretion rate is 
indicated by a dotted line and time is taken in units of sec- 
onds. The shock wave moves up and down over the spatial 
pre-shock zone, but it maintains the Mach number of ~ 2 
throughout the shock oscillation. 

Although in model B with the highest A the shock is 
formed firstly around r ~ 100i? g and moves gradually up- 
stream, it finally disappears and a steady flow is attained. 
In model C with A smaller than model A, the flow features 
are similar to model A. But the shock oscillates between 
40-80 R g and the mass-outflow rate varies by a factor of 4. 
Model D differs from the other models because the initial 
flow of model D is based on the ID adiabatic flow with sta- 
tionary standing shocks. Fig. 4 shows the time variations of 
L, Mout, Mi n and R s for model D. First, the shock is formed 
steadily at r/R g ~ 20 and it begins to move between 15 
< r/Rg < 30. The luminosity is modulated by several per 
cent but the mass-outflow rate varies by a factor of 7 al- 
though it is two orders of magnitude smaller than the input 
accretion rate. The Mach number of the shock maintains it- 
self at about ~ 4 and is stronger than those in models A 
and C, because the shock is formed closer to the black hole. 

We may wonder that the flow of model D shows non- 
stationary behaviours, because the initial flow of model D 
has been given from the adiabatic solution with the station- 
ary standing shocks. So we examined this model under an 
adiabatic 2D condition and reconfirmed that the stationary 
shock is formed constantly at r ~ 16R g which is slightly 
smaller than the theoretical values of R s — 17 and 19 ob- 
tained from the ID adiabatic solution. In model D, unlike 
the adiabatic case, we notice that the cooling and heating 
of the gas considerably influence the formation and position 
of the shock. 

Fig. 5 shows the time variations of the total luminosity 



L by the synchrotron and free- free emissions, M ou t, Mm and 
R s for model F. The figure also shows the irregularly oscil- 
lating shock phenomena similar to model D but the modu- 
lation amplitude of the luminosity is much larger than that 
of model D and the maximum luminosity is as high as five 
times of the averaged luminosity 1.6 x 10 36 erg s _1 which is 
comparable to the observed luminosity of Sgr A*. In the two- 
temperature models of E and F, the synchrotron emission 
is one orders of magnitude larger than the free-free emis- 
sion and the energy transfer rate q lc from the ion to electron 
by Coulomb collisions is negligibly small compared with the 
synchrotron emission. 

Table 3 shows the averaged values of L, M[ n / Mi nput and 
Mout /Mi n put, the shock position R B on the equatorial plane, 
the maximum luminosity L max and the maximum mass- 
outflow rate (M ou t)max, where Mi nput is the input accretion 
rate of 2.5 x 10 20 g s _1 . The luminosities of the models A-D 
are ~ 5 x 10 38 erg s" 1 independently of A. The radiative 
efficiency r\ (= L/Mc 2 ) of the accreting matter into the ra- 
diation is ~ 10 -3 in models A-E but is far low, ~ 10~ 5 , in 
model F. The mass-outflow rate increases with increasing A. 
In model B with the highest angular momentum, most of 
the input accreting matter is ejected as the wind, without 
arriving at the inner edge, and in model D with the lowest A, 
only four per cent of the input matter is ejected as the wind. 
Conversely, the mass-inflow rate at the inner edge decreases 
as A increases. 



5.2 Structure of the flow and of the shock 

Fig. 6 shows the profiles of density p (g cm~ 3 ), gas tem- 
perature T (K), radial velocity v, its Mach number on the 
equatorial plane and accreting mass flux M(r) (g s^ 1 ) ver- 
sus radius at t = 1.7 x 10 5 s for model A. The accreting mass 
flux is defined as the mass flux of the flow with negative ra- 
dial velocities above the equatorial plane. Since the flow is 
convectively unstable in the inner region of r ^ 30R g , M(r) 
is not correctly defined at r < 30i? g according to the above 
definition. The non-stationary shock is found at a radius of 
r ~ 90i? g . The temperature 2.5 x 10 9 K at the outer bound- 
ary increases to 5.0 x 10 9 K at the pre-shock radius and is 
enhanced to 1.3 x 10 10 K across the shock. The density is 
also enhanced by a factor of three across the shock. In Fig. 6, 
a single shock structure is found but we find two shock fea- 
tures at some phases. The accreting matter falls inward with 
a constant accretion rate until it attains the shock position 
and then it decreases (see the scale on the right side of the 
figure) after the gas passes through the shock. This follows 
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Table 3. Averaged values of L, M; n /M; nput and M out /Mi nput , shock position R a , maximum luminosity L max and maximum mass- 
outflow rate (Mout)max, where Mi nput is the input mass accretion rate. 



Model 


L (erg s 1 ) 




Mout /M input 




imax / L 


(M ou t) max /M out 


A 


5.4X10 38 


0.32 


0.68 


80-160 


1.04 


~ 3 


B 


5.3xl0 38 


1.0xl0~ 4 


0.99 


no shock 






C 


5.1xl0 38 


0.72 


0.28 


40-80 


1.08 


~ 4 


D 


5.1xl0 38 


0.96 


0.04 


15-30 


1.05 


~ 7 


E 


1.3xl0 38 


0.30 


0.70 


70-130 


~ 2 


~ 2 


F 


1.6xl0 36 


0.97 


0.03 


7-45 


- 5 


~ 7 



from the fact that the outflow begins in the post-shock re- 
gion after matter crosses the shock. Fig. 7 shows the profiles 
of p (g cm -3 ), T (K), v, its Mach number and M(r) (g s" 1 ) 
at t — 2.0 x 10 5 s for model D. The shock on the equatorial 
plane is found at r ~ 17J? g . In this figure, we also find that 
the accreting mass flux M(r) becomes a little small in the 
post-shock region although it is almost constant until the 
gas reaches the shock. This is due to the effect of the mass- 
outflow in the post-shock region. However, the effect on the 
accreting mass flux is smaller than in model A, because the 
mass-outflow rate in model D is far smaller than in model 
A. 

Fig. 8 shows the profiles of density p (g cm -3 ), ion tem- 
perature Ti (K) , electron temperature T c (K) , radial velocity 
v, Mach number of the velocity on the equatorial plane and 
accreting mass flux M(r) (g s" 1 ) for model F at t = 1.3 x 10 6 
s, at which the luminosity takes a maximum value 3.2 x 10 38 
erg s _1 . The profiles in the pre-shock region are almost same 
as the initial one and the ratio T e /Ti is not so different from 
its initial ratio 8(= 1/9) even in the post-shock region. This 
means that the steady values of the electron and ion temper- 
atures and accordingly the intrinsic total emission of models 
E and F may not be still reached by the time integration 
over t ~ 2 x 10 s s. To obtain exact steady distribution of 
T e and Ti, the time-independent differential equations for a 
stationary state should be solved instead of the method of 
the time-dependent calculation. In spite of this fact, from 
the result of the two-temperature model with the input ac- 
cretion rate of 4.0 x 10 -6 Mq yr -1 , we can expect the total 
emission of ~ 10 36 erg s _1 compatible with the observed one 
of Sgr A* if the electron temperatures of ~ 10 9 - 10 10 K are 
obtained in the inner region. 

Figs 9 and 10 show the contours of density p (g cm -3 ) 
and temperature T (K) with velocity vectors at t = 1.7 x 10 5 
s in model A. The shock wave rises obliquely at r ~ 90i? g 
above the equatorial plane and extends upward to ~ 50i? g . 
The velocity vectors in Fig. 10 are amplified by five times, 
compared with Fig. 9. The initial accreting matter at the 
outer boundary supersonically falls down towards the gravi- 
tating centre and is decelerated at the shock front, enhancing 
the density and the temperature and becomes subsonic. The 
post-shock region with high densities and high temperatures 
results in a strong outward pressure- gradient force also along 
the Z-axis and a part of the accreting matter deviates from 
the disc flow and escapes as the wind flow. The wind velocity 
becomes 0.01-0. 05c at the outer boundary except for a nar- 
row funnel region along the Z-axis, where the funnel wall 
is characterized by the vanishing of the effective potential 
and the flow velocity within the funnel region is relativistic 



at ~ 0.2c. The remaining accreting matter is swallowed into 
the central black hole. In model A, two thirds of the input 
accreting matter is expelled from the accretion flow. The 
mass-ouflow begins after the matter undergoes the shock 
compression. The circulating flow is also found at r $C 30i? g , 
that is, the accretion flow is convectively unstable in the 
post-shock region. 

5.3 Properties of the time-variability of the flow 

The most remarkable feature of the low angular momentum 
flow considered here is the existence of the non-stationary 
and oscillating shocks found in all models except for model 
B. This shock phenomenon results in the irregular modula- 
tion of the flow variable. . 

The luminosity varies synchronously with the mass- 
outflow rate without any time lag. In spite of the large 
modulation of the shock position in models A, C, and D, 
the modulation amplitude of the luminosity is as small as 
a few to ten per cent , while in models E and F it varies 
by a factor of 2-5. The outflow consists of a persistent out- 
flow and a large episodic one with irregular intervals of an 
hour to ten days and the episodic outflow rate is higher by 
a factor of 2-7 than the averaged mass-outflow rate. 

Sgr A* possesses a quiescent state and a flare state. The 
flares of Sgr A* have been detected in multiple wavebands 
from radio, sub-millimetre, IR to X-ray l|Yusef-Zadeh et aD 
I2009l . l201lf ). The amplitudes of the variability at radio, IR 
and X-ray are by factors of below 1/2, 1-5 and 45 or even 
higher, respectively, a few times per day. The rapid vari- 
ability in the radio band on time-scales of a few se conds 
to hours has been observed l|Yusef-Zadeh et ah! 1201 ll ). Al- 
though a quasi-periodicity of ^ 20 min u tes in the IR flares 
dGenzel et al.ll2003l : iMever et al]|2006al lbl; iTrippe et af]|2007l ; 
lEckart et alJ 120081) and radi o QPOs of 16, 22, 31 and 56 
minutes l)Mivoshi et alJlioill ) have been reported, they are 
the s ubject of a controversy l|Mever et alj|2008l ; |Po et al~l 
2009). These flare phenomena are interesting in relation to 
the time variations of the luminosity and the mass-outflow 
rate in the present models. We consider that the intermit- 
tently large modulations of the outflow due to the oscillatory 
shock may be related to the flares of Sgr A*. Unfortunately, 
the time resolution in our calculations is at most as small 
as a few minutes because of the limited time step in the 
calculations. Therefore, the short time variability at a scale 
of ~ 10 minutes observed in Sgr A* can not be compared 
with the result of our models. Figs 11 and 12 show the power 
spectrum density of the total luminosity L and mass-outflow 
rate M out and of the disc luminosity L d and shock radius 
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Figure 3. Time variations of to- 
tal luminosity L, mass-outflow 
rate Mout, mass-inflow rate M; n 
and shock position R B on the 
equatorial plane for model A, 
where input accretion rate is in- 
dicated by a dotted line and time 
is taken in units of seconds. 
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Figure 4. Same as Fig. 3 but for 
model D. The luminosity is mod- 
ulated by several per cent but 
the mass-outflow rate varies by a 
factor of 7 although it is two or- 
ders of magnitude smaller than 
the input accretion rate. 



R s on the equatorial plane, respectively, for model A, where 
the disc luminosity is defined at the surface of the inward 
flow with only negative radial velocities near the equatorial 
plane. Here, we can not find a conspicuous frequency peak 
in these power spectra partly due to insufficient statistical 
power but they are suggestive of semi-regular variabilities 
on time-scales of hours to days. 

The origin of the oscillating shock appearing in the 
present models is entirely attributed to the range of param- 
eters of the specific energy and the specific angular momen- 



Figure 5. Time evolutions of total emission L by the synchrotron 
and free-free emissions, mass-outflow rate M ou t, mass-inflow rate 
M; n and shock position R B on the equatorial plane in model F. 



turn used here. Generally, the shock formation and its posi- 
tion are very sensitive to the upstream flow parameters. The 
initial flow variables of the models, except for models D and 
F with A = 1.35, are based on the ID adiabatic flow in which 
parameters A never produce a stationary standing shock but 
are located just nearby in the range of parameters responsi- 
ble for the stationary shock. Furthermore, we know that the 
shock becomes unstable w hen we take account of the energy 
loss and gain of the gas jMolteni et al.lll996l ; lOkuda et al.l 
l2004h . 



6 COMPARISON WITH THE 

ADVECTION-DOMINATED ACCRETION 
FLOW 

The present low angular mome ntum flow models ar e com- 
pared with the ADAF mode l dNaravan fc Yil Il994 Il995l; 
Naravan fc McClintockl 120081 ; lYuan. Quataert fc Naravanl 
2003h which is referred to as a radiatively inefficient ac- 
cretion flow model of Sgr A. The region considered here 
is limited to an inner region of r ^ 200i? g . The ADAF has 
generally the following properties: (1) the accretion flow be- 
comes geometrically thick, (2) the pressure support in the ra- 
dial direction is considerable, so the angular velocity is sub- 
Keplerian, (3) the radial velocity of the gas is much larger 
than that of the standard disc model and is barely less than 
Keplerian velocity in the innermost region, and (4) the large 
radial velocity and the large scale height cause the gas den- 
sity to be very low, and so the cooling time is very long com- 
pared with the dynamical time and the matter is optically 
thin. These properties are also realized in the present low an- 
gular momentum flows. The ADAF solution is gas-pressure 
dominated and the sound velocity is comparable to the Kep- 
lerian velocity, so the gas temperature becomes nearly virial. 
Gas at such a high temperature radiates copiously at small 
radii where the temperature can approach 10 12 K. Thus, in 
the ADAF model, a two-temperature plasma of the ion tem- 
perature Ti and the electron temperature T c is assumed to 
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r/R g 

Figure 6. Profiles of density p (g 
cm ), temperature T (K), ra- 
dial velocity v, Mach number of 
the radial velocity on the equa- 
torial plane, and accreting mass 
flux M(r) versus radius at t = 
1.7 X 10 5 s for model A. 
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Figure 7. Same as Fig. 6 but 
at t = 2.0 X 10 5 s for model 
D. The parameters A and e in 
this model correspond to ones re- 
sponsible for stationary standing 
shocks under the adiabatic con- 
dition. 

exist at the small radii. Typical ADAF models have th e two 
temperatures scaling as l|Naravan fc McClintockll2008r ) 

7] ~ \Q 12 /r K, T c ~ Min( 2], 10 9 " 11 K). (17) 

In our models A-D, the two-temperature model is not con- 
sidered but a high temperature of ~ 2 x 10 11 K comparable 



Figure 8. Profiles of density p (g cm -3 ), ion temperature T\ (K), 
electron temperature T c (K), radial velocity v and Mach number 
of the radial velocity on the equatorial plane versus radius at 
t = 1.3 X 10 6 s for model F. 



to the above ion temperature is obtained in the innermost 
region. The temperature profile is roughly T oc r" 1 in all 
models, in agreement with the ADAF. In the ADAF with 
the constant accretion rate, a self-similar flow with the den- 
sity profile of p(r) oc r~ s gives s = 1.5, while our models 
show s = 1.4-1.8. The radiative efficiency of the accretion 
flow in our model F is as small as ~ 10~ 5 in the usual 
ADAFs. 

Thus, as far as the inner region of the accretion flow is 
concerned, the low angular momentum flow model has many 
properties in common with those of the ADAFs. The distinc- 
tion between the ADAFs and the present low angular mo- 
mentum flow is that the former describes a hot, nearly spher- 
ical, accretion flow with viscosity, while the latter shows the 
same accretion flow with lower angular momentum but with- 
out viscosity. The remarkable difference between our models 
and the ADAFs is the existence of the oscillating shock phe- 
nomena in the low angular momentum model, except for 
model B with a relatively larger angular momentum. The 
shock position R s is irregularly oscillatory in the inner re- 
gion of r ^ 180-Rg. Another distinctive point of our models 
is that the mass-outflow is considerable, depending on the 
specific angular momentum of the model, as mentioned in 
Section 5.1. 

Fig. 13 shows the energy spectra calculated from model 
F with L = 3.2 x 10 3e erg s" 1 at t = 1.3 X 10 6 s (dashed 
line) and with L = 1.2 X 10 36 erg s" 1 at t = 7.1 x 
10 5 s (solid line). The sub-milimetre and X-ray bumps in 
these spectra are originated from the synchrotron cooling 
by thermal electrons and the bremsstrahlu ng cooling, re- 
specti vely, whose formulas are p r escrib ed in iNaravan fc Yil 
l|l995h and IStepnev fc Guilbertl (|l98St ). The intensities in 
the sub-milimetre bumps of models E and F are a little 
higher and the peak positions lie in a longer wavelength com- 
pared with the observational spectra at the quiescent state 
of Sgr A*. To the contrary, the intensities in the X-ray band 
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Figure 9. Contours of density p 
(g cm -3 ) and velocity vectors at 
t = 1.7 X 10 5 s for model A. The 
velocity of light is indicated by 
the upper arrow. The shock wave 
extending obliquely is formed at 
r I R g ~ 85 up to a height of z ~ 
50R g . 




r/Rg 



Figure 10. Same as Fig. 9 ex- 
cept for the temperature, where 
the velocity vectors are ampli- 
fied by five times, compared with 
Fig. 9. 

(2-8 keV) are one orders of magnitude lower than the ob- 
served one. The radiatively ineffic ient accretion flow models 
ijYuan. Quataert fc Naravarj2003r i. such as the ADAFs, take 
into account the complicated physics including synchrotron 
and inverse Compton emission by thermal and non-thermal 
electrons under the two-temperature assumption and de- 
scribe the energy spectra very fitted to the observed ones 




v 

Figure 11. Power spectrum 
density of the total luminosity L 
and the mass-outflow rate M ut 
in model A. 
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Figure 12. Same as Fig. 11 but 
for disc luminosity and the 
shock position R s on the equa- 
torial plane. 

of Sgr A* over a wide range of energy bands from the radio 
to X-ray. In this respect, the present our model is too sim- 
ple and insufficient to reproduce the spectra of Sgr A* but 
it is useful to examine the basic scenario of the low angular 
momentum flow model. 



7 SUMMARY AND DISCUSSION 

Based on the analysis of the angular momentum of the ac- 
cretio n flow around Sgr A* by Moscibrodzk aTTjas fc Czernvl 
(2006), we examined a low angular momentum flow model 
for Sgr A* using two-dimensional hydrodynamical calcula- 
tions. The models with the specific angular momenta A of 
1.35, 1.55, 1.68, and 2.16 (in the usual nondimensional units) 
and input accretion rates of Minput = 4.0 x 10 -6 was con- 
sidered. We summarize the results as follows. 
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Figure 13. Energy spectra vL v erg s of model F with L = 
3.2 X 10 36 erg s^ 1 at t = 1.3 X 10 6 s (dashed line) and with 
L = 1.2 X 10 36 erg s" 1 t = 7.1 X 10 5 s (solid line) 



(1) The accretion flow is very hot, optically thin, geo- 
metrically thick, and convectively unstable in the inner re- 
gion of the accretion flow, as is also found in the ADAFs. 

(2) The flow is highly advected and, if the syncrotron 
cooling and a two-temperature model are taken into account, 
the radiative efficiency of the accreting matter into radia- 
tion is very low, ~ 10 _5 -10 -3 , and the input accretion rate 
Minput = 4.0 x 10~ 6 Mq yr _1 results in the total luminosity 
~ 10 36 erg s _1 which is comparable to the observed lumi- 
nosity of Sgr A*. 

(3) An irregularly oscillating shock is formed in the in- 
ner region of a few tens to a hundred and sixty Schwarzschild 
radii, except for the model with the highest A = 2.1 6. Due to 
the oscillating shock, the luminosity and the mass-outflow 
rate are modulated by several per cent to a factor of 5 and 
a factor of 2-7, respectively, on time-scales of an hour to ten 
days. 

(4) The mass-outflow rate increases with increasing A 
and it ranges over a few per cent, one-thirds, two-thirds and 
99 per cent, of the input accreting matter, corresponding to 
the cases of A = f .35, f .55, f .68 and 2.f6, respectively. 

(5) The oscillaing shock is necessarily triggered if the 
specific angular momentum A and the specific energy e be- 
long to or are located just nearby in the range of parameters 
responsible for a stationary shock in the rotating inviscid 
and adiabatic accretion flow. 

(6) The time variability may be related to the flare phe- 
nomena of Sgr A*. 

Since the specific angular momenta used in the models 
are reasonably estimated, the time variability of Sgr A* may 
be attributed naturally to the oscillating shock in the low 
angular momentum model considered here. Further evidence 
of the low angular momentum of the accretion flow around 
Sgr A* would confirm the validity of the models. If the spe- 
cific angular mommentum of the accretion flow around Sgr 
A* is actually as small as A = 1.35, the time variability itself 
of the total emission like model F may explain the quiescent 
and the flare states of Sgr A*. On the other hand, in the case 
of higher angular momentum i.e. A = 1.68, the time mod- 



ulation of the total emission L is small to explain the flare 
phenomena of Sgr A*. However, under this circumstance , 
we expect that the large episodic wind caused by the non- 
stationary shock may induce the flare phenomena of Sgr A*. 
The flare phenomena of Sgr A* have been see n to be polar- 



ized in the radio, sub-mil li metre and IR b ands (Gcnz el et al 



20031; iMever et a l. 2006a b]; iTrippe et af]|2007l ; 



Eckart et al 



2006). This is a strong evidence for a magnetic field around 
Sgr A*. The episodic wind in our model may be accelerated 
as a relativistic jet by some mechanism such as magnetic 
reconnection process in a current s eet under the ma gnetic 
field away from the accretion flow |Yuan et al.l l2009f). The 
jet will collide with the ambient matter and excite the mat- 
ter as a hot emitter, especially at times of the episodic wind. 
The emission from the hot matter may be observed as the 
flares in X-ray. Furthermore, the jet expands and cools in 
the far distant interstellar matter and results in the syn- 
chrotron emission due to the electrons in the magnetic field. 
This may be observed as the radio flares with a time lag 
from the X-ray flares. The existence of the magnetic field 
plays an important role in the accretion flow. However, the 
magnetic field was not taken into account simultaneously 
with the hydrodynamics in our models. Thus, the current 
low angular momentum flow model is simple and can not 
reproduce the observed spectra of Sgr A*, as is successfully 
fitted over a wide range of the wavelength by ADAFs. For 
a spectral fitting, we have to treat the physics of the syn- 
chrotron and inverse Compton emission by the thermal and 
non-thermal electrons under the two-temperature model. A 
further magneto-hydrodynamical examination of the flow 
with a low angular momentum model including the above 
physics should be pursued in the future. 
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